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Abstract. Gravitational collapse is analyzed in the Brane-World by arguing that regularity 
of five-dimensional geodesics require that stars on the brane have an atmosphere. For the 
simple case of a spherically symmetric cloud of non-dissipating dust, conditions are found for 
which the collapsing star evaporates and approaches the Hawking behavior as the (apparent) 
horizon is being formed. The effective energy of the star vanishes at a finite radius and the 
star afterwards re-expands and “anti-evaporates”. Israel junction conditions across the brane 
(holographically related to the matter trace anomaly) and the projection of the Weyl tensor 
on the brane (holographically interpreted as the quantum back-reaction on the brane metric) 
contribute to the total energy as, respectively, an “anti-evaporation” and an “evaporation” term. 


1. Introduction 

Black holes (BHs) are unstable in four (and higher) dimensions because of the Hawking effect |I], 
which is deeply linked to the trace anomaly of radiation fields j2J. In the Randall-Sundrum (RS) 
Brane-World (BW) models [3], a collapsing homogeneous star likewise requires a non-static 
exterior [1]. Further, forcing a static exterior induces a trace anomaly of the same form as that 
of semiclassical BHs, although with opposite sign, which suggested that BH solutions of the bulk 
equations correspond to quantum corrected (semiclassical) BHs on the brane (5J f|, in the spirit 
of the holographic principle El and AdS/CFT conjecture JH|. 

The junction conditions j2], which preserve the regularity of (five-dimensional) geodesics, 
cannot allow a step-like discontinuity (e.g. across the star surface) and a (5-like discontinuity 
(e.g. across the brane) in the stress tensor at the same location, hence discontinuities in the stress 
tensor of brane stars are not mathematically permitted. This can be physically understood by 
considering that the brane thickness (of the order of the AdS length t ~ A -1 / 2 , A being the brane 
tension) and the star’s atmosphere cannot be both negligibly thin, and it is indeed natural to 
assume that the latter is much larger than l. In Ref. m, we employed effective four-dimensional 
(hydrodynamical) equations m in order to study a “corrected” Oppenheimer-Snyder (OS) 
model m in which the star is divided into three regions (see Fig. 0]): a homogeneous “core” 
with most of the energy; a “transition region” of fast density decrease which connects with a 
“tail”, where the energy density slowly vanishes. The tail and transition region together form 
the “BW atmosphere”, which disappears in the limit A —> oo. 


In general, effective equations cannot determine the brane metric uniquely unless one also 
knows the bulk geometry. However, for negligible dissipation and asymptotically flat brane, 
the evolution of the system is uniquely determined by the dynamics of the homogenous core, 
which corresponds to an exact five-dimensional solution m and reproduces the trace anomaly of 
quantum field theory J2]. Our main results are that the total energy of the system is conserved 
and that the collapsing star “evaporates” until the core experiences a “rebound” in the high 
energy regime (when its energy density is comparable with X), after which the whole system 
“anti-evaporates”. This behaviour cannot be related to GR perturbatively (in e = po/X , where 
po is the initial core density), but it seems in agreement with the uncertainty principle of quantum 
mechanics ESJ. Moreover, we can find a range of parameters for which the minimum radius of 
the collapsing core is larger than t (which sets the scale of Quantum Gravity in the BW), thus 
supporting the holographic interpretation of the model. 


2. Spherically symmetric collapsing dust 

Effective four-dimensional Einstein equations for dust and vanishing brane cosmological constant 
can be written as [ll] = 8 tt T®® = 8 7r ^ p cS u^ u u + , where u M is the unit 

four-velocity of dust, h^ the space-like metric that projects orthogonally to iG, 
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with p the (“bare”) energy density, and U and come from the Weyl tensor projected on the 
brane 1 . Bianchi identities supplied by the junction conditions then produce both local (LCE) 
and non-local conservation equations (NLCE). For negligible dissipation, one can take a Tolman 
brane geometry, ds 2 = —dT 2 + ( R'y d?’ 2 + i ? 2 dfl 2 where R = R(r, r), and the LCEs imply 
conservation of the “bare” mass function 

rn p( r ) = f Q P(E x ) <9x (R 3 (t,x)) dx . (2) 

The system of NLCEs is in general not closed, since we do not have an equation for 
However, for sufficiently large R, the knowledge of Td^ v in an extended spatial region together 
with the asymptotic flatness and the continuity of the Weyl tensor make that system closed. 
Another important result which follows from the LCEs and NLCEs is that, if the brane metric 
is asymptotically flat, the anisotropic stress / 0 whenever p' 0. 

The (r, r) Einstein equation yields the equation of motion 
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(3) 


where the (in general time-dependent) “effective” mass 

M(r,r) = ~^- f Q p eS (r,x)d x (ii 3 (r,a:)) dx , (4) 

replaces the bare mass of General Relativity (GR). The latter is always well defined and dust 
shells in GR move along geodesics to reach the central singularity (R = 0) at increasing proper 
times (Tolman model EH) or at the same proper time (OS model [T2|). In the BW, the effective 
mass M drives the collapse, but it diverges for R —* 0 and this would make the whole four¬ 
dimensional brane singular. In order to avoid this (physically unlikely) case, one has to include a 
sufficiently negative contribution to the mass from the projected Weyl tensor which will generate 
an Hawking flux near the forming horizon and make the effective mass evaporate completely. 

1 We set the momentum density Q^ = 0, see Ref. m for more details. 
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Figure 1. Density profile. 
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Figure 2. Qualitative behavior of V for 
p > p c (solid line) and p < p c (dashed line). 
For p = p c (dotted line), the peak of V equals 
the shells energy E = 0. 


2.1. Core 

The bulk solution which corresponds to the OS core of the star is regular in five dimensions El 
and, since p' = 0, the system of relevant equations is now closed with II = 0 and 
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where p is a constant. The physical radius R can be written in the factorized form 


R(r,r) 
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in which Mg is the total bare mass of the OS core, and the effective mass © is given by 


M (t, r ) 


M s 


(-) 3 + ^- 
VO/ 32 7T Pq 



(2M S ) 2 /r\ 
3 R 3 {r 0 J 


R 


(7) 


where the second term represents a BW correction. The above mass would diverge for R —* 0, 
the harmless (at least when covered by an horizon) central singularity of GR, and make the 
whole space-time singular (this also occurs for a more general Tolman core uni). One must 
therefore have p positive and large enough so that each shell will bounce back after reaching 
a minimum radius R m [ n where the corresponding effective mass vanishes. The Weyl tensor, 
holographically interpreted as the quantum back-reaction on the brane metric, contributes the 
“evaporation” term proportional to p which dominates at low energies; the BW correction to 
the matter stress tensor, holographically interpreted as the matter quantum trace anomaly 0, 
yields the “anti-evaporating” term proportional to Mg which increases with the energy. 

Upon inserting the effective mass © into Eq. ©, one obtains 
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which shows that the core remains “rigid” through the bounce and it is thus sufficient to consider 
the evolution of its surface at r = r^. There exists a critical value p c = (32 n po Mg/3e)^^ 3 such 
that the bounce occurs only for p > p c (for p = p c the two turning points coincide, see Fig.©. In 
Fig. ©we display a typical trajectory of Rq(t) = R(t, rg), along with the corresponding effective 






















Figure 3. Typical evolutions of the core 
radius Rq(t)/Ro{0) (solid line) and effective 
mass Mo(t)/Mo(0) (dashed line) for p > p c . 
Time is arbitrary. 



Figure 4. Energy flux versus the time t = 
(r — r§ s )/T (solid line) compared to the 
Hawking flux (dotted line) for e = 10~ 4 , Mg = 
po = 1, T = 10 and p = 5000 > 4824 = p c . 


mass Mo(r) = M(r, ro), for p > /i c . After the rebound, the system reverses its evolution, 
however, in a more realistic model, the collisionless description of dust should be relaxed and 
dissipation from both the core and the atmosphere is expected to make the process irreversible. 
Since Rq < 0 and Ro(0) 3> 2 Mg (there is no initial horizon), for p > p c 
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and the evaporation sets out at the beginning of the collapse and goes on at least until the 
radius bounces back at R m i n rsj l (Mg/£) 1/3 A 1//2 . Since i? m ; n is the shortest length in 

our system, the holographic description is expected to hold for astrophysical sources for which 
Mg 3> A -1 / 2 (it also follows that p c 1). 

The horizon trajectory r = ?’h(t) is defined as the value of r at which the horizon is formed 
at the time r, that is R(t,th(t)) = 2 M(r, rn(r)). The time derivative of the horizon effective 
mass, Mh(t) = M(r,rn(T)), 

= M„ + m , (10) 

then contains two contributions: one, from the intrinsic time dependence of M at constant r, 
is a first order effect in e which vanishes in GR; the second one, due to the (possibly) variable 
number of shells included within the horizon, depends on the detailed form of the atmosphere. 

Since in the core, the velocity \R(t, r)| increases monotonically in r at fixed r, an (apparent) 
horizon forms at the boundary r = ro at r = Ty S = T — | Mg, where T fixes the time scale of 
the collapse. From Eq. ©, we then get the (instantaneous) Hawking flux [1; 
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precisely at the time r = Tjj S . 


2.2. Transition region 

In the GR model p = 0 for r > ro- The density must therefore decrease rapidly from p = O(e 0 ) 
at r = ro to p = O(e) for r = r s . Moreover, since the transition is a BW effect, we can take 
f s — ro = 0(e) as well as U = 0(e). Combining these results, we obtain, to first order in e, 


M(r, r) ~ M 0 (r) , 


(12) 



















for ro < r < r s , and the Hawking flux will remain negative (and substantially unaffected) 
throughout the border of the transition region r = r s . 


2.3. Tail 

As in the transition region, p = 0(e) for r s < r, and p' p/X = 0(e 2 ), so that bulk gravitons 
are decoupled from brane matter. The Weyl contribution is however of the same order, 
rnu(T, r; r s ) = ^ J)T U ( R 3 )' da; = O(e), and satisfies 
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Since R(r,r) > Ro(t) for r > ro, (fl3l) implies that mu cannot remain zero in the tail. 

In order to proceed, we now assume that: 

(i) Hindoo R(t, r) = oo, and 

(ii) the effective mass be finite at spatial infinity , Hindoo M(r,r) < oo, Vr > 0. 

Since Mq(t) always remains finite if p > p c and the tail bare mass is small by construction, this 
implies that lim^oo mu(j, r; r s ) < oo, Vr > 0. Asymptotic flatness ensures that we can take 
the limit r —> oo (equivalent to R —> oo at fixed time) in Eq. (11 HI) and finally obtain 

lim muir, r;r s ) + Mo( r ) = Am M(r, r) = 0 , Vr > 0 . (14) 

r—>oo r— kx) 


We have thus shown that if the total effective mass at spatial infinity is finite at the initial time 
t = 0, it will always remain constant (for a bouncing core evolution with p > p c ), so that the 
total effective mass of the collapsing dust star is actually conserved. 

It is particularly interesting to consider the case for which there is initially no energy stored in 
U. We then see from Eq. m that mu > 0 during the collapse and, after the bounce, we expect 
that mu will also evolve backwards so as to ensure the condition ®. The energy flux along 
the horizon trajectory r = rn(r) can also be obtained numerically for this case and is plotted 
in Fig. for some values of the parameters. We can see that it grows less with respect to that 
predicted by Hawking and the conclusion is that, although the evaporation sets out according 
to Hawking’s law, the back-reaction on the brane metric subsequently reduces the emission. 


3. Luminosity 

A distant observer experiences an impinging flux of energy during the collapse, 
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which shows the same dependence on the mass Mh as the semiclassical expression when the 
horizon is first forming, and subsequently decreases to zero (before it becomes negative). 
However, since this happens after the apparent horizon begins to form, a distant observer might 
have to wait an infinite amount of time to measure a vanishing flux. 

Upon equating the BW result m at the time when the horizon starts to form to the 
semiclassical luminosity calculated in the Schwarzschild background for an astrophysical object, 
we obtain a bound for the AdS length m , 

10~ 32 mm <C l <S 10“ 9 mm , (16) 

which is three orders of magnitude better than the best constraint found in m considering the 
time scale of primordial BH evaporation. 









4. Trace anomaly 

Strictly speaking, there is no trace anomaly in our approach, since we have included the back- 
reaction of the effective matter on the brane metric. However, in order to compare with known 
results without the back-reaction, we can define the trace anomaly 7Z as the sum of the Ricci 
scalar and the trace of the bare stress tensor. At the OS boundary, r = ro, we then have 

_ 9 Me 

1Z = -h. 
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which is the quantum Ricci anomaly of Ref. 0 with the correct sign. It is then clear that the 
sign mismatch found in Ref. [3j was due to the choice of a non-smooth energy density and that 
the atmosphere describes how the excess energy stored in the OS boundary is released. 

5. Conclusions 

Inspired by the conjecture that classical BHs in the BW may reproduce the semiclassical behavior 
of four-dimensional BHs, we have studied the gravitational collapse of a spherical star of dust 
in the RS scenario. Regularity of the bulk geometry requires continuity of the matter stress 
tensor on the brane and leads to a loss of mass from the boundary of the star. We found 
that the system of effective BW equations is closed to our level of approximation and leads to 
the collapsing dust star emitting a flux of energy which, at relatively low energies, approaches 
the Hawking behavior when the (apparent) horizon is being formed. However, no real space- 
time singularity forms since the star effective mass vanishes at finite star radius, thus leaving a 
remnant which re-expands by absorbing back the emitted radiation. In a more realistic case, one 
however expects that dissipation cannot be neglected and the process then becomes irreversible. 

Let us finally point out that all the above features were obtained for BHs formed by 
gravitational collapse, excluding therefore primordial BHs. 
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